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This  paper  deals  with  the  design  and  analysis  of  an  adaptive  array 


processor  for  use  in  passive  detection  of  directional  stochastic  signals 


such  as  are  found  in  sonar.  The  processor  consists  of  a set  of  tapped  delay 


line  filters,  one  for  each  array  element.  The  algorithm  used  for  adjusting 


the  tap  gains  is  a modification  of  the  stochastic  approximation  method  of 


Robbins  and  Monro,  and  it  utilizes  knowledge  of  the  signal  autocorrelation 


function  and  spatial  direction  of  the  target.  It  is  shown  that  the  final 


form  approached  by  the  processor  is  that  of  a space-time  filter  optimized 


in  the  direction  corresponding  to  assumed  target  location 


The  system  performance  is  analyzed  by  considering  a noise  field 
consisting  of  a spatially  isotropic  component  and  a single  directional 


component  referred  to  as  an  interference.  It  is  shown  that  although  the 


system  initially  cannot  discriminate  between  the  target  and  the  inter 


ference,  it  eventually  acts  to  eliminate  the  effect  of  the  interference 


almost  completely.  It  is  also  shown  how  the  useful  signal  to  noise  ratio 


increases  during  adaptation 
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I.  Problem  Statement  and  Objectives 

The  problem  considered  in  this  paper  is  the  passive  detection  of  a 
noise-like  signal  waveform  generated  by  a source  located  in  a known  direction 
from  the  receiver.  Typical  applications  of  this  general  problem  can  be 
found  in  sonar  detection,  seismic  detection  and  radio  communications.  The 
sonar  application  is  the  one  that  primarily  motivates  this  study,  and 
examples  will  be  taken  from  the  sonar  area.  In  order  to  take  advantage  of 
the  known  directivity  of  the  target  signal,  a directional  receiver  in  the 
form  of  an  array  is  employed  to  distinguish  signal  from  noise.  In  the 
sonar  application,  the  receiver  consists  of  an  array  of  hydrophones, 
together  with  an  appropriate  processor.  Generally  speaking,  the  processor 
consists  of  individual  filters  on  each  sensor  output,  a summer,  a post- 
summation filter,  a square-law  device  and  an  averaging  filter.  The  output 
of  the  averaging  filter  is  used  to  indicate  the  presence  of  a target 
signal. 

In  the  absence  of  a target  signal,  the  averaging  filter  output  is 
the  result  of  noise  waveforms  picked  up  by  the  array  elements.  The  noise 
is  partly  far-field  noise  and  partly  locally  generated.  The  far-field 
noise  is  often  assumed  to  be  directionally  isotropic;  however,  there  may 
also  be  directional  noise  sources.  These  directional  noise  sources  are 
referred  to  as  interference  sources  while  the  directionally  isotropic 
component  is  referred  to  as  ambient  noise.  In  the  absence  of  interference 
noise,  detection  of  a target  signal  can  be  based  simply  on  the  presence  of 
a directional  component  in  the  received  signal.  However,  if  interference 
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sources  can  be  expected  to  be  present  in  the  noise  field,  then  it  is  lutt  se am 

necessary  to  define  the  target  signal  in  some  way  to  distinguish  it  from  toll*rrn.l 
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This  paper  is  concerned  with  developing  a system  for  processing  the 
outputs  of  a passive  array  of  hydrophones  under  the  following  assumptions: 

(1)  Target,  interferences  and  ambient  noise  are  assumed  to  be  Gaussian 
random  processes. 

(2)  The  sum  of  interferences,  ambient  noise,  and  local  noise,  are  regarded 
as  the  effective  noise,  which  is  assumed  to  be  statistically  indepen- 

J 

dent  of  the  target  signal. 

(3)  The  target-signal  component  s^(t)  observed  at  the  output  of  the  i*"*1 
hydrophone  is  a linear  time-invariant  transformation  of  d(t),  the 
target-signal  component  observed  at  the  output  of  an  ideal  isotropic 

. hydrophone  located  at  the  origin  of  the  coordinates.  The  target 
direction  is  known,  together  with  its  autocorrelation  function  (but 
not  necessarily  its  power  level). 

(4)  The  statistics  of  the  noise  field  are  completely  unknown.  Interferences 

may  be  present,  but  this  is  not  known.  If  they  are  present,  their 

/ 

directions  are  unknown. 

(5)  The  wavefronts  of  target  and  interferences  are  regarded  as  plane  over 
the  dimensions  of  the  receiving  array. 

The  proposed  system  consists  of  an  adaptive  linear  multichannel 
tapped-delay-line  filter  and  algorithms  for  iterative  adjustment  of  the 

filter  coefficients 

on  the  tapped-delay  lines.  The  adjustment  algorithm  is  based  on  the 
method  of  stochastic  approximation,  and  it  utilizes  the  knowledge  about 
target  direction  and  autocorrelation  function  that  is  assumed  to  be 
available.  Thus  a space-time  filter  optimum  in  a predetermined  direction, 
and  designed  to  suppress  interference  signals  from  other  directions  is 
produced.  By  varying  the  azimuth  for  which  the  filter  is  optimized  the 
system  produces  a bearing  response  pattern  which  can  be  examined  by  an 


2. 


1 


i 


li 


Mh 


operator  to  determine  whether  a target  is  present.  In  many  practical 
situations,  narrow  peaks  might  be  interpreted  as  targets,  while  broader 
peaks  might  be  classified  as  interference. 


3. 


II.  Array  Systems 

The  receiving  array  is  assumed  to  consist  of  K omnidirectional  receiving 
elements  (sensors).  Each  element  has  a total  output 

*i(t)  = s±(t)  + ni(t)  (1) 

consisting  of  a signal  component  s^(t)  plus  an  undesired  noise  component 
ni(t). 

Consider  a general  array  configuration  consisting  of  individual 
filter  H^io)  on  each  sensor  output  x^t),  i = 1,2,...K,  a post-summation 
filter  G(w),  a square- law  device,  and  an  averaging  filter  h (to).  The 

3V 

output  of  hav(uj)  is  used  to  detect  the  presence  of  a target  signal  in  noisy 
environment.  See  Figure  1. 

Let  <t^(w)  be  the  signal  power  spectral  density,  a_(w)  be  the  steering 

vector  with  elements  expCju)-^),  i = 1,...K,  and  let  $^(0))  be  the  spectral 

1 2 3 A 

matrix  of  the  input.  It  has  been  shown  by  many  authors  ’ ’ ’ that  the 
optimum  filters  are  given  by 
H («)  = [®*  (w)]"1  a*(u) 

P ^ (2) 

G(oj)  = 4>  ^(oj) 

for  Gaussian  processes  and  small  input  signal-to-noise  ratios.  The  whole 
system  forms  a likelihood  ratio  detector  and  at  the  same  time  the  output 
signal-to-noise  ratio  is  maximized. 
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Ill . Adaptive  Delay  Line  Filters 

In  the  proposed  system  the  individual  filters  H^(u>)  are  constructed 
in  the  form  of  tapped-delay-line  filters  consisting  of  a tapped  delay  line, 
adjustable  gains  whose  input  signals  are  the  signals  at  the  delay-line 
taps,  a summer  to  add  the  weighted  signals,  and  machinery  to  automatically 
adjust  the  gains.  The  impulse  response  of  such  a discrete  system  is 
completely  controlled  by  the  weight  settings. 

th 

Referring  to  Figure  2,  let  c^  be  the  weight  at  the  k 1 tap  on  the  i 
filter,  and  t^  the  time  delay  at  the  k**1  tap.  Let  the  delayed  signal  be 

\k  * sik + vik  ‘ *i(t  ‘ 

where  £ ^ is  the  signal  component  and  v ^ the  noise  component.  Define 

- = ^C10C11***C1MC20*“C2M"‘CK0"‘CKM1 

(3) 

- = ^n10nll*’,nlMn20’*,n2M",nK0*,*riKM^ 

The  mean-squared  error  (mse)  between  the  desired  cr  target  signal 
d(t)  and  the  summer  output  z(t),  is  minimized  by  setting  the  weights  to 

w = (R  + R r1  R,  = R-1  R,  (4) 

-op  — £ — v — d£  — n — d£ 

where  R^  and  Ry  are  respectively  the  delayed  signal  and  noise  correlation 

matrices  of  dimension  K(M  + 1)  x K(M  + 1) . R. _ is  the  correlation  function 

-dL 

vector  between  the  desired  signal  and  the  various  delayed  signals. 

The  stochastic  approximation  algorithm  used  to  adjust  the  filter  weights 
takes  the  form 

-j+1  = ~j  + 2yj^d£  “ 2yjzjnj 

Y 1 

In  Eq.(5),  the  y's  satisfy  y = — - — with  — < a<_  1 and  R^  is  completely 

j ^ 

specified  once  correlation  function  and  angular  direction  of  the  signal 
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are  known.  A block  diagram  of  the  adaptive  mechanism  Is  shown  In  Figure  3. 
The  convergence  properties  of  algorithm  (5)  have  been  studied  in  a 

t; 

previous  paper  and  more  extensively  in  Ref.  6.  It  has  been  shown  that 

algorithm  (5)  converges  in  mean  square  and  in  probability  as  long  as  the 

second  order  statistics  of  the  input  processes  are  bounded  and  if  of 

Eq.(4)  is  a constant  vector.  The  rate  of  convergence  depends  on  the  input 

statistics,  various  system  parameters,  and  training  environment.  The 

inverse 

mean-squared  error  is  found  to  decrease  approximately  as  the/first  power 
of  the  adaptation  time.  The  rate  of  convergence  is  essentially  independent 
of  the  number  of  weights  to  be  adjusted  as  the  algorithm  allows  simultaneous 
adjustments.  The  size  of  the  error,  however,  depends  on  the  total  number 
of  taps  and  the  starting  point. 

In  order  to  determine  the  perf romance  of  the  adaptive  receiver,  the 
mean  values  of  the  weights  during  the  adaptation  period  are  derived  here 
for  reference. 


Since  the  summer  output  in  Figure  2 is 
T 

z.  = w.  ' n (6) 

j “3  ~3 


r 1 


and  a new  delayed  input  vector 

n'  ■ P n. 
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Since  R^  = ^n”op»  as  seen  from  Eq.(4),  we  can  transform  Eq.(9)  into 


- (1  - 2y.A)  E[w'  ] + 2y  A w* 

J J “ j j op 


(12) 


or 


E(^+l>  - Kp  ■ (1  - 2'j  - Hip)  (13) 

Now  consider  a particular  component  of  w',  and  for  clarity  no 
subscript  or  superscript  denoting  the  component  is  used.  Then  we  obtain 
a difference  equation  for 


- “iP  ■ (1  - 2y><"j  - »;p> 

whose  solution  is 


w 


j+1 


= (wi  ~ woP)  n (1  ~ 2y<a)  + *1 


k-l 


op 


(14) 


(15) 


If  we  choose 


j 2 (j+1) 

[This  means  that  Yj  is  the  jth  element  of  a diagonal  matrix  having 
K(M  + 1)  elements],  and  note  that 


(16) 


J 


J 


n (i  - 2y  x)  - n (i  - t-~ ) 

k=l  * k=l  k+1 


1 

j+1 


we  have 


Wj+1  “j+1  ^ + j+1  wo 


op 


(17) 


(18) 


Although  the  choice  of  Yj  given  by  (18)  results  in  the  simplest  expression 
f°r  wj+i>  ot^er  choices  give  similar  results.  See  Ref.  6. 
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Some  Remarks  on  the  Operations  of  the  Proposed  System 


(a)  Choice  of  the  Initial  Weights 

Although  the  adjustable  weights  can  be  set  to  any  values  at  the 
beginning  of  the  adaptive  process,  it  is  desirable  to  set  them  not  too 
far  from  their  optimum  by  using  whatever  information  is  available  concerning 
the  statistics  of  the  noise  field.  The  formula  for  calculating  the  optimum 
gains  can  be  utilized  to  start  the  initial  computation  with  inaccurate  noise 
statistics.  This  kind  of  choice  will  shorten  the  adaptation  period  and 
thus  reduce  the  cost  of  operation.  In  cases  where  absolutely  no  such 
information  is  known,  the  gains  associated  with  the  input  delays 
T_^(i  = 1,2,...K)  are  set  to  1 and  the  rest  to  zero.  This  results  in  a 
conventional  power  detector.  As  the  adaptation  proceeds,  the  whole  system 
will  gradually  be  transformed  into  an  optimum  one. 

(b)  Problem  of  Signal  Suppression 

In  adaptive  detection  systems  such  as  those  proposed  by  Glaser^, 

g 

or  Jackowatz  , where  adaptation  depends  on  signal  information  produced  by 
the  processor,  there  is  a critical  input  signal-to-noise  ratio  below  which 
the  system  rejects  the  signal.  This  does  not  take  place  in  our  system  since 
the  signal  information  utilized  by  the  algorithm  is  that  supplied  by  R,  , and 
this  is  independent  of  input  SNR.  The  adaptation  process  basically  pro- 
duces a space-time  filter  that  is  optimized  in  a given  direction  regardless 
of  whether  a signal  is  actually  present  in  this  direction  or  not. 

(c)  Unknown  Signal  Power 

Although  it  may  be  reasonable  to  assume  knowledge  of  the  shape 
of  the  signal  spectrum,  it  would  generally  be  unrealistic  to  assume  anything 


_ to  be  known  about  the  power  level.  It  is  shown  in  Appendix  D of  Ref.  6 that 
if  the  assumed  signal  power  differs  from  the  actual  value  by  a multi- 
plicative constant  then  the  filter  weights  converge  to  a value  multiplied 
by  that  same  constant.  Since  a constant  multiplier  applied  to  a filter 
leaves  the  filter  unchanged  (except  if  its  output  is  compared  to  some  fixed 
threshold)  lack  of  knowledge  of  the  signal  power  level  causes  no  real 
change  in  the  basic  system. 


IV.  Performance  of  the  Adaptive  Receiver 

Our  major  concern  is  to  deal  with  the  passive  detection  of  a sonar 
target  in  the  presence  of  ambient  noise  as  well  as  interferences  from  other 
targets. 

In  order  to  illustrate  the  essential  procedures,  we  make  the  following 
simplifying  assumptions: 

a)  The  array  is  linear  with  elements  spaced  a constant  distance  d apart. 

b)  The  input  spectra  are  identical  in  shape  over  a frequency  range 
(0,1*0;  above  the  frequency  the  signal  power  is  zero. 

c)  The  array  dements  are  sufficiently  far  apart  so  that  the  ambient 
noise  is  statistically  independent  from  element  to  element. 

d)  There  is  only  a single  interference  source. 

The  mean  value  of  the  detector  output  in  the  presence  of  the  target 
signal  is  easily  found  to  be 


S+N  2 n 


— / |G  | 2 HT$  H*du> 

tr  1 1 — vv — 


where  the  explicit  dependence  of  the  integrand  on  to  is  omitted  for 
simplicity.  A similar  expression  is  obtained  for  <y>^,  the  mean  value  of 
the  output  in  the  absence  of  target  signal,  by  replacing  by  j>  . 
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Hence  the  d.c.  change  of  the  output  due  to  the  presence  of  the  target  is 


v>s+N  - <*>,,-  IT  ' l°l2  »\/d“ 


For  small  input  SNR  the  output  variance  can  be  shown  to  be  approximately 

o 2 = — ^ — / | G | 4 { HT  <t>  H*}2dw  (21) 

y it  T — nn 

■ av  -00 

where  is  the  averaging  time.  The  output  SNR  is  then  defined  by 

y , /o  . 
dc  y 

Let  0^,0^.,  and  0 be  respectively  the  target  and  interference  bearing 

and  the  array  steering  angle.  These  angles  are  measured  relative  to  a 

line  perpendicular  to  the  array  axis. 

Let  t = ~ sin  0_  (22) 

o c T 

P = 7 sin  0 (23) 

o c I 

^ = - sin  0 (24) 

o c 

where  d is  the  hydrophone  spacing  and  c is  the  sound  velocity.  0 is  the 

independent  variable  in  the  directivity  pattern  yg+^(6)*  Then  the  signal 

and  interference  delays  at  the  i1"*1  array  element  are 

t = (K  - i) t (25) 

i o 

Pt  = (K  - i)pQ  (26) 

and  the  steering  delay  is 

t = (K  - i)^  (27) 

1 o 

Let  the  interference  delay  vector  be  given  by 
T jwp,  jtop_  jwp 

b^u)  - [e  e e k]  (28) 

Then  because  of  assumptions  (c)  and  (d)  at  the  beginning  of  this  section 

<t>  (w)  = d>  (u>)U  + <}>_  (u)  b b*T  (29) 

- nn  n — i 


where  U is  the  identity  matrix  and  where  <j>  (w)  and  <}>  (w)  are  the  ambient 

n I 

and  interference  noise  spectral  densities  respectively.  Also 

$ (w)  = <t>  (w)  + <t>  (w)  = <j>,(w)  a a T + $ (to)  (30) 

Note  now  that  if  optimum  expressions  for  H are  substituted  in  Eqs . (19), 

(20),  and  (21)  the  scalar  spectral  functions  (w) , (to) , and  <t>^(u))  appear 

only  as  ratios  <p,(u)/<p  (w)  <(>,.  (w) /<{>  (u)  etc.  Thus  the  assumption  that  the 

an  in 

spectral  shapes  of  signal,  noise,  and  interference  are  identical  permits 

us  to  replace  these  ratios  by  constant  ratios  S/N,  I/N,  etc.  Also,  since 

the  signal  power  is  assumed  to  vanish  for  frequencies  beyond  , the 

integrals  in  (19) , (20)  and  (21)  are  evaluated  between  the  limits  -w  to  u . 

o o 

The  effect  of  assumption  (b)  is  therefore  the  same  as  if  signal,  ambient 

noise,  and  interference  spectra  were  assumed  to  have  flat  spectra  of 

height  S/2,  N/2  and  1/2  respectively  in  the  frequency  range  [-wo,wo],  and 

to  be  zero  otherwise.  This  assumption  is  the  one  that  is  actually  used  below. 

til 

Suppose  now  that  the  transfer  function  of  the  i filter  is 
M 

H (w)  = E c e (31) 

X k=0 


Then  Eqs.  (19),  (20),  and  (21)  become: 
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<y= 


S+N 


S 

4tt 


K 

E 


K 

E 


M 

E 


w 


1 cikch*  ; dw  ^X  X 

i=l  h=l  k=0  H=0  1 -w  XiXh 

o 


ju(A 


h£ 


Aik) 


(32) 


_ 2 K K M M o jw(x  - x ) jw(A  - A ) 

y.  - fr  E E E E c.,  c,  f / du>  e 1 n e M lk  (33) 
d’C*  8*  i=l  h=l  k=0  £=0  ik  U -a, 


and 


2 KKKKMMMM 
o = E E E E E E E E 
y i=li'=l  h=l  h'=lk=0  k’=0«,=0  «.'=0 


CikCi 'k ' Ch£Ch ' 1 ' 


0) 


4irT 


o 

/ dm  9 


Jo<V  + Vf  ‘ 4ik  - Vk'> 


av  -a) 


ninh  ni'nh’ 


(34) 
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In  these  expressions  0 (or  0 ) is  the  ij  element  of  the  input  (or 

Vh  ninh 

noise  spectral  matrix  and  Eq.(34)  is  valid  for  the  case  of  small  signal- 
to-noise  ratios. 

(a)  Initial  Behavior 

Assuming  the  worst  case  where  absolutely  no  information  about 
the  noise  field  is  known,  the  gains  associated  with  the  input  delayed  by 
x^(i  = 1,2,... K)  are  set  to  1 and  the  rest  to  zero  so  that  a square-law 
detector  is  used  at  the  start.  Hence  the  output  of  each  array  element  is 
delayed  to  provide  maximum  response  in  the  signal  direction,  i.e., 

-i(w)  = N -*  (35) 

The  weights  and  spacings  are  simply 

Cik  = 6ik/N 

(36) 

Aik  = Vik 

where  6..  is  the  Kroneker  delta, 
ik 


Substituting  Eq.(36)  into  Eq.  (33)  gives  the  d.c.  change  of  the 


output 


1%  £ 

* N2 


where  we  have  assumed  that  the  integral  of  e % 0 for  A ^ 0.  Similarly, 


the  output  variance  is 

2 r 

o K io  . . _ , K-l  sin  to  ip 

a 2 = (i)2  1 + I [2  + A j.  o__o_ 

y T it  V 1 N U K . . u ip 

J=1  av  _ i=l  o o 

t2  . K-l  sin  or  ip 

+ ? (1  + 1 * 

N i=l  o o 

„ K-l  K=1  sin  to  (i-h)p  sin  w (i+h)p 
. ±_  T r r o O O 0 

if 2 w (i-h)p  or  (i+h)p 

K i-i  n=l  oo  o o 


(K-i) ] 


] (K-i) (K-h) } 


Dividing  Eq.  (37)  by  the  square  root  of  Eq.(38)  gives  the  output  signal- 
to-noise  ratio,  which  becomes  approximately  [10] 


12. 


- T co  1/2  „ T T2  „ , -1/2 

SNRl4^  f KU+  2 N + ~2  (f  K+  3K")} 

N 


(39) 


if  the  maximum  frequency  processed  is  very  high  such  that  w p » 1 and 

o o 


the  terms  associated  with  pQ  make  negligible  contributions  except  for  i = h. 


For  most  cases  of  practical  interest,  the  number  of  hydrophones  in 
an  array  is  large,  K >>  1,  so  that  for  ambient-noise-dominated  environment 


SN^  « K(|)  when  (|)2  » | K 


(40) 


and  for  interference-dominated  environment 


SNR1  “ K1/2(|)  when  (|)2  <<  -|  K 


(41) 


The  average  output  of  the  squarer,  y,  yields  the  so-called  directivity 
pattern,  which  may  be  obtained  by  varying  the  electrical  time  delays. 


In  the  signal  direction  x = x and  for  to  p >>  1 Eq.(32)  becomes: 

o o o o * 


>i<8-  V :s?KlI1  + l + Kfi 


(42) 


In  the  interference  direction  x = p , hence 

o o 


y'l  ( 0 " 0T>  §[(!  + £)  Hi] 


'1'  ” “I 
and  in  any  other  direction 


S, 

N' 


¥ 


(43) 


yi(e)  ^K|d  + i + l> 


(44) 


These  expressions  indicate  quite  clearly  that  there  are  peaks  in  the 


directivity  pattern  in  the  signal  and  interference  directions  (6  = Q^and 


0 = 0^)  and  that  the  system  does  not  discriminate  in  any  way  against 

the  interference. 


- - ■ - ■■  > ■ ^ - - 
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(b)  Final  Behavior 
i.  Optimum  Gains 

Since  the  adaptive  processor  converges  to  the  optimum 

detector,  the  final  values  of  the  gains  are  given  by  Eq.(4).  If  the  input 

signal- to-noise  ratio  is  small  R £ R ; also  it  can  be  shown6  that  if 

-T1  “ M 

all  the  spectra  are  similar  R = K,R  where  K,  > 1 is  a constant.  Hence 

~ M 1“"V  1 — 

we  may  use 

W„  = r;1  R^  (45) 

This  equation  is  not  suitable  for  obtaining  analytical  expressions  for 
the  optimum  gains,  and  therefore  we  approximate  the  gains  by  a simple 
Fourier  integral  expression: 


cik  = YT  /Hi(“)e  d“  (46) 

.00 

where  H^(w)  is  the  i1"^  optimum  filter  given  in  Eq.(2)  with  4^  replaced 

by  $ 

— nn 

Then,  using  Eq. (29)  in  (2)  we  find  that  the  final  expression  for  the 
entire  filter  is  given  by 

H »(»)  - O'1*/72  & ■ 4-  '2  - iTTSTf1  i <"> 


The  i row  of  Ha>(w)  is 

* K * 
b . £ b,  a. 

V“>  ■ T t»i 5:5 1 

K + N/I 


so  that  the  impulse  response  is 

K jco(t-pi+pk-Tk) 

h l ^ : r jw(t"ti)  k=ie 

hi(t)  =ir~r  1 [e ; 

K + N/I 


■]  dw  (49) 


It  is  clear  that  the  impulse  response  of  the  individual  filters  can  be 
completely  specified  by  setting  the  gains  according  to 


C = — (6 — T-) 

ik  N ^ ik  K + N/l' 


at  delays 
A.,  = 


pi  - °k  + \ 


Note  that  since  a practical  delay  line  can  supply  only  discrete  delays  it 
is  generally  not  possible  to  obtain  the  exact  values  required  by  (51).  It 
is  assumed  that  the  delay  lines  have  a sufficient  number  of  taps  so  that 
the  performance  degradation  produced  by  this  is  negligible, 
il.  Signal-to-Noise  Ratio 

If  we  substitute  Eq. (50)  into  Eqs.(33)  and  (34),  and  follow  the 
same  steps  leading  to  Eq.(39),  we  obtain  the  signal-to-noise  ratio  of  the 
detector  output: 

cmr  - I 1/2  i K(K-l)  N/l 

“ " 2 C 7T  ' N K + N/I 


, K-l  sin  to  ip 

i * i (K-i)  — 

K(K-1  N/I)2  i=l  “oipo 


„ K-l  K-l 

~2 ? 1 1 (K-D  (K-h) 

K (K-1+N/I)  i=l  h=l 


sin  (o  (i-h)p  sin  w (i+h)p 

o o o o 

to  (i-h)p  ..... 

oo  to  (i+h)p 

o o 


>>  1,  then 

i 

i T to  1/2 

It  av  °) 

Z 7T 

|(K-1) 

. T hi  1/2 

2 1 ^ °) 

|(K-D 

1+N/KI 


3 K(K-l+N/l) 2 


Equation  (52)  or  (53)  gives  the  asymptotic  performance  of  the 
adaptive  array  processor.  Since  the  training  period  is  finite,  the  actual 


fr 
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signal-to-noise  ratio  is  lower  than  that  given  by  Eq.(52)  or  (53).  As 

would  be  expected  Eq. (52)  is  just  the  output  signal-to-noise  ratio  of  an 

9 

optimal  (likelihood  ratio)  detector  first  investigated  by  Schultheiss 

and  then  by  Tuteur^  from  a simpler  formulation,  but  under  the  same 

assumption  of  similar  input  spectra  over  (0,io  ). 

o 

iii.  Directivity  Pattern 

Although  the  directivity  pattern  can  be  obtained  by  using  Eq.  (50) 
in  (32),  the  fact  that  (50)  represents  essentially  the  optimum  filter  makes 
it  possible  to  use  Eq.(47)  in  (19)  to  put  the  result  into  a more  compact 
matrix  notation: 


a)  . *T 

— IS  ° A *t  — — 

y°o(e)  = YiT  N f dm  a [U  - K+N/J 


b b 


*T 


ny  + |bb*T  + f aa^ltU-^a  (54) 


-ID 


where  it  is  to  be  understood  that  the  a/s  and  b/s  are  functions  of  to. 

Expansion  of  the  integrand  results  in  twelve  terms  typified  by  the  term 
~*X  *T  *T  - -2 

a b.  b b^  ji  (K  + N/I)  . The  vector  product  that  remains  after  re- 
*T 

moving  b b = K is 


K 
Z 

i=l  h=l 


K -p  +e  -t  ) 

E e 


We  consider  the  integral  of  the  summand  to  be  approximately  zero  except 

for  i = h,  in  which  case  it  is  2to  . Thus  the  contribution  of  this  term 

o 

2 

is  2K  to^.  The  other  eleven  terms  are  evaluated  similarly  except  that  in 

^ *x  *T  ~ 

combination  such  as  a a a a or  a bb  a one  must  consider  also  whether 

A 

the  independent  variable  tq  is  near  t or  near  pq  respectively.  Thus,  if 
* ~*T  *X~  2 *• 

t =i,a  aa  a = K , if  t differs  from  T the  contribution  of 
o o ’ o o 

'*T  *T- 

a a a a is  only  K. 

By  utilizing  considerations  of  this  sort  we  obtain  the  following 
approximate  expressions  for  the  directivity  pattern: 
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a)  For  angles  for  away  from  target  or  interference 


y 2.  5.  / 1/  i _i_  / N/I  s2  A . KN/I  A S ,r„  , , , N/I  ,2 


ya>(6)  % f-  f (K-l  +(K+N/I 


) + 


(K+N/I) 


2 + N [K_1  + (K+N/l 


y) 


(55) 


In  the  signal  direction  (0  = 0^) 


1'*,<e  ■ eT’  ^ !r  f <K_1  + <K+?7r')2  + 


KN/I  S f 2K  2K  - 1 

(K+N/!) 2 N " N/I+K  (K+N/!)2 


]} 


and  in  the  interference  direction  (0 


V 


- /fl  _ n ^ '0°  S N/I  .2  K2N/I 
yco(6  0I)  ^ Tt  N ^k^k+n/i  + 


_ . S ( N/I  2. 
2 + N K(K+N/I)  1 


(57) 


(K+N/I) 

Note  that  if  the  noise  is  ambient-dominated  N/I  >>  K,  and  in  this  case  these 
three  expressions  are  essentially  the  same  as  Eqs . (42) , (43) , and  (44).  This 
is  to  be  expected  since  in  the  absence  of  interference,  and  if  all  spectra 
have  the  same  shape  the  simple  filter  of  Eq.(35)  is  optimum.  On  the  other  hand, 
if  the  interference  noise  is  large  N/I  <<  K,  the  three  equations  above  become 
respectively: 

c c 

(55a) 


y-«»)  *^f(K-D(f+  i) 


y«>(0  = eT)  % ^f(K-i)(K  |+i) 


<v,  o S v r /N/I \ 2 , S ^ ^ (N/I)' 


700(0  = 0j)  ^1TWK[(  K } § + « + 


(56a) 

(57a) 


Thus  the  optimum  filter  tends  to  suppress  the  interference  noise  in  all 
directions,  and  for  large  I/N  it  produces  a dip  in  the  response  pattern  in  the 
interference  direction. 

(d)  Adaptive  Behavior 

Duringthe  adaptation  period  the  transfer  function  of  the  i^  filter 
is  given  by  Eq.(31).  The  c's  appearing  in  this  equation  are  represented  by 
the  vector  w of  Eq.(3),  which  changes  with  time  according  to  Eq.(5).  The 
variation  of  the  elements  of  w during  the  training  period  therefore  determines 


(56) 
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the  adaptive  behavior  of  the  processor.  Since  in  Eq.(5)  and  are 
random,  the  w^  are  also  random.  Thus  we  consider  only  the  expected  values 
of  w . According  to  Eq.(18)  the  expected  values  of  w^  is  related  to  the  initial 
and  final  value  by  the  general  formula 

EI«j+1l  “ + (58) 

where  w,  is  the  initial  value  of  wand  w the  final  value,  and  where  the  p. 

—1  — —co  3 

and  q^  are  functions  of  j that  depend  on  the  choice  of  the  weighting  sequence 
Yj  of  the  stochastic  approximation  algorithm.  It  is  clear  that  p^  = 1, 
qL  = 0,  p„  - 0,  qK  - 1. 

Combining  Eqs.(31)  and  (58)  we  can  write 

«j+1  - + qjHooCw)  (59) 

where  H^(u)  and  Hco(<i))  denote  respectively  the  initial  and  final  forms  of 
the  transfer  function  vectors.  For  the  particular  example  considered  here, 

H^(u))  is  given  by  Eq.(35)  and  H x(u)  by  Eq.  (2). 

Basically,  we  are  required  to  evaluate  the  following  three  integrals 
(*) 


2 7! 

/ 

0 

1 

IT  T 

av 

1 

o ^ 
3 

cj+i  • rv  ; lGl2(Pj»i + «ij»(-)Tixx<pJ«1  + q-ji^  (62> 

where  in  Eq. (61)  the  H's  use  the  steering  vector  a rather  than  the  signal 
delays  a. 

Note  that  in  general 

T *2T*2T* 

(p.H.  + q H 4>(p  H + q,H  ) = p.4H,  ♦ H.  + q.  H 4>  H 

rj— 1 j °°  — 1 j — 1 1 j 


+ 2pjqjH1  ♦ H 


i 
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The  first  two  terms  of  Eq.(63)  are  already  available  from  the 
sections  on  the  initial  and  final  behaviors.  Thus,  if  we  set 


previous  two 


3 (i) 

Aj+1  = i=l  Aj+1 


K2u>  2 

where  y by  Eq.(37).  Similarly 

j-1  v N 


(64) 


j-1 


(2)  2 1 ,i„i2  „T  - * 


A.v:;  = q.“  ■—  I\g\*  HA  4>  H*  dm  = q . 2y  , 
j+1  nj  2 TT  _oo  1 ' -«  — ss— 00  Mj  ^d.c. 


(65) 


A 


(3) 

j+1 


CO 

2p.q,  -7T-  /|G|  VT  * H*  dw 
11  2 u 1 1 —1  -SS  00 

J J oo 


K2oj 


2p,q 


£ /Sv2  K-l  + N/I 

'AT  ' 


j \j  TT  'N'  K + N/I 


(66) 


Similar  expressions  can  be  obtained  for  BJ , , and  C,.,;  B.t,  requires  the 

j+1  j+1  J+1 

evaluation  of  six  terms,  the  evaluation  of  three  terms.  For  details  see 

[6].  Then  the  output  signal-to-noise  ratio  is 


SNR 


= _J±1 
3+1  7b  )1/2 

(Bj+i) 


(67) 


and  the  directivity  pattern  is 


yj+l(e)  Cj+1  (68) 

Eqs.  (65)  and  (66)  have  been  computed  by  means  of  a digital  computer,  using 
the  optimum  y given  by  Eq.(16),  Figure  4 shows  the  variation  of  output  signal- 
to-noise  ratio.  Figure  5 shows  the  variation  in  the  directivity  pattern 
for  target  and  interference  well  separated  in  bearing.  The  way  in  which  the 
interference  is  suppressed  during  the  adaptation  is  clearly  evident  in  Figure  5. 
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